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We study one-dimensional quantum systems near the classical limit described by the Korteweg- 
de Vries (KdV) equation. The excitations near this limit are the well-known solitons and phonons. 
The classical description breaks down at long wavelengths, where quantum effects become dominant. 
Focusing on the spectra of the elementary excitations, we describe analytically the entire classical- 
to-quantum crossover. We show that the ultimate quantum fate of the classical KdV excitations is 
to become fermionic particles and holes. We discuss in detail two exactly solvable models exhibiting 
such crossover, the Lieb-Liniger model of bosons with weak contact repulsion and the quantum 
Toda model. We argue that the results obtained for these models are universally applicable to all 
quantum one-dimensional systems with a well-defined classical limit described by the KdV equation. 
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I. INTRODUCTION 


The Korteweg-de Vries (KdV) equation 1 describes 
propagation of waves in a medium with competing dis¬ 
persion and nonlinearity and is ubiquitous in the physics 
of classical nonlinear systems.® In addition to the usual 
periodic waves, the KdV equation supports solitons, i.e., 
localized disturbances with particle-like properties. Pro¬ 
tected by the integrability of the KdV equation, the soli- 
tons move and scatter off each other without distortion. 

Real systems are never integrable. Nevertheless, the 
KdV equation often provides an adequate effective de¬ 
scription of the excitations in the long-wavelength limit, 
whereas deviations from the integrability have a signif¬ 
icant impact on the lifetime and stability of solitons at 
shorter wavelengths.® Another limitation on the applica¬ 
bility of the KdV equation, which is the primary focus 
of this paper, arises in one-dimensional quantum systems 
with a well-defined classical limit. In these systems quan¬ 
tum effects inevitably become dominant at sufficiently 
long wavelengths, leading to the breakdown®!!!] 0 f the 
classical description. 

Because the emergence of quantum behavior is not 
associated with broken integrability, it is natural to 
approach the problem of describing the classical-to- 
quantum cr ossov er from the perspective of exactly solv¬ 
able models. 1 ®® Fortunately, solvable models exhibiting 
such crossover are available. In this paper, we consider 
two well-known examples. The first one is the Lieb- 
Liniger model of bosons with weak contact repulsion. 1 USI 
In the classical limit, its dynamics can b e de scribed by 
the mean-field Gross-Pitaevskii equation, 15 16 which has 
the form of the nonlinear Schrodinger equation. 2 -! The 
equation is integrabl^® and reduces to the KdV equa¬ 
tion in the long-wavelength limit P' T he se cond model we 
consider is the quantum Toda model.® i^SUlts classical 
counterpart, the Toda modelp®is also integrable and the 
corresponding equation of motion in the continuum limit 
coincided^ with the KdV equation. 

It is well known that at high momenta the exact spec¬ 


tra of elementary excitations of the Lieb-Liniger and the 
quantum Toda models m atch those deduced from their 
classical analogs.! 12 ! 18 ! 22 ! In this paper we focus on low 
momenta, and evaluate the spectra near the classical-to- 
quantum crossover. Some of the results presented below 
have been reported in Refs. |SI9j . 

The rest of the paper is organized as follows: In Sec. [Tl] 
we summarize some well-known facts about the KdV 
equation and discuss its applicability for the description 
of quantum interacting systems near the classical limit. 
In Sec. Ill we introduce the Lieb-Liniger and the quan¬ 
tum Toda models and discuss their relation to the KdV 
equation. Spectra of the elementary excitations are eval¬ 
uated in Sec. EY1 The results are discussed in Sec. m 
Technical details are relegated to the Appendixes. 


II. KDV EQUATION IN QUANTUM SYSTEMS 


The KdV equation in the standard form read;P 

drf + 6/0 £ / + dlf = 0, (2.1) 

where r and £ are the dimensionless time and coordinate, 
respectively. The equation can be viewecP as a canon¬ 
ical equation of motion d T f = {f,TL} generated by the 
Hamiltonian 

U = jdi 

and the Poisson bracket 

{/(O./(O} = -0e*(£-O- (2.3) 


/ 3 - 2^/) 2 


( 2 . 2 ) 


The KdV equation supports an infinite number of poly¬ 
nomial integrals of motion,® all in involution with each 
other with respect to the Poisson bracket (2.3). Among 
these integrals of motion are the dimensionless momen¬ 
tum 


V = 



(2.4) 
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and the Hamiltonian (2.21 itself; the latter defines the 
dimensionless e nergy . Substituting the single-soliton 
solutiorP of Eq. (12.1 ), 


/o(£,t) = 


2A 2 


cosh 2 [_4(£ — 4*4 2 r)] ’ 


(2.5) 


into Eqs. (2.2) and (2.41 and excluding the parameter A. 
we obtain the relation 


Ho 


1 ^ 37^0 ^ 5/3 


( 2 . 6 ) 


between the corresponding values of the dimensionless 
energy and momentum. 

In this paper, we study one-dimensional interacting 
quantum systems with a well-defined classical limit de¬ 
scribed by the KdV equation. As we will demonstrate in 
Sec. |III[ in these systems the low-energy right- and left- 
moving excitations decouple from each other. Focusing 
on this regime and on, say, the right-moving excitations, 
we introduce a bosonic field >I> obeying the commutation 
relation 


[4>(;r), $(y)] = i7T sgn(cc — y) (2.7) 

and write the Hamiltonian of the right-movers as a sum 
of two contributions, 

H = vP+ H KdYl (2.8) 

where 

P= ^ f L dx:{d x ^>) 2 -. (2.9) 

4ttJ 0 

is the momentum operator and 


PfCdV — 


h 2 


127nn* 


dx : 


(d x $) 3 


( 2 . 10 ) 


represents the leading perturbation. In Eqs. (2.9) 


and (2.10), L is the size of the system, and the colons 


denote the normal ordering with respect to the bosonic 
vacuum. The parameters v, to* and a* have the units 
of velocity, mass and length, respectively; their origin 
and physical meaning will be elucidated below. We are 
interested in the excitations of the Hamiltonian (2.8) 


with wavelengths of order a*. Provided that the length 
scale a* is small compared with L , spectra of such excita¬ 
tions are not sensitive to the precise form of the boundary 
conditions imposed on the field $(x), which we take to 
be periodic, 


4?(x) = 4>(a; + L). 


( 2 . 11 ) 


Comparison of Eqs. (2.7), (2.9), and (2.10) with 


Eqs. (2.3), (2.4), and (2.2), respectively, shows that, 


apart from rescaling, d x <& is merely the quantized ver¬ 
sion of /. The classical model defined by Eqs. (2.2) and 


(2.3) and its quantum versiorP^I defined by Eqs. (2.7) 
and ( 2.10| represent the so-called first Hamiltonian struc¬ 
ture of the KdV equation. 2 The quantum KdV model 
is believed to be integrable, with the first few integrals 
of motion having the same forrrPS! as the correspond¬ 
ing classical expressions. These models should be distin¬ 
guished from those corresponding to the second Hamil¬ 
tonian structure, in which both the Hamiltonian and 


the Poisson bracket have forms different from Eqs. (2.2) 
and ( |2.3| .^4l Whereas quantum KdV models of the latter 
type have received much attentiorP® 2 ' due to their rele¬ 
vance to conformal field theorypSl the former arise natu¬ 
rally as the low-energy description of one-dimensional in¬ 
teracting quantum systems, see, e.g., Refs. j5l9linl29l30j . 

Because the operators P and Pi<dv commute, the ex¬ 


citation spectrum of the Hamiltonian (2.8) has the form 


e(p) = vp + e Kdv {p), 


( 2 . 12 ) 


where e, p, and eicdv are eigenvalues of the operators H, 
P, and -ffKdV; respectively. In order to find the nonlin¬ 
ear in p contribution SKdv(p) in Eq. ( 2. 12| ), we consider 
the Heisenberg equation of motion <9 t $ = MjT^-ffKdVj d>], 
which after a change of variables 


X — a*^, t — t*7”, 

takes the form 


1 * - 


3 to* a 2 


(2.13) 


d r 4> + -:(<%$) 2 : +9|4> = 0. 


(2.14) 


We now differentiate Eq. (2.14) with respect to £ and 
treat the field 


F =\ <%$■ 


(2.15) 


as a classical variable. Neglecting quantum fluctuations 
of F amounts to replaci ng (: F 2 :) with ( F ) 2 and leads 
to the KdV equation © for the expectation value 
/ = (F). Solution of this equation satisfying the condi¬ 
tion fo^ a *dt; /(£,r) =0 [see Eqs. (2.11) and (2.15)] can 
be written as 


/(£,t) = -Af 0 ^- + fo{£,r), 


(2.16) 


where A/"o = / d£ = 2(3Po/2) 1 ^ 3 is the dimension¬ 

less mass of the soliton^ and / 0 is given by Eq. (2.51. 


We now take the expectation values of P and Pkuv 


see Eqs. (2.9) and (2.10)], neglect quantum fluctuations, 


and compare the resulting expressions with Eqs. (2.2) 
and (2.4). In the limit a*/L —► 0, this yields 


P= -Poj ^Kdv — ——'Ho- (2-17) 

TTCL^ 71’Cjjj 

Upon introducing the momentum and energy scales 


P* = 


3h 

2a* 


27 H* 

£* =- 

8 t. 


2 TO* 


(2.18) 
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we write the KdV contribution to the spectrum as 

^Kd v(p) = £*e{p/p*). (2.19) 

For the soliton excitation the dimensionless function e(s) 
in Eq. (2.19) is given by 


, X 3 (2w\ 

^soliton $) — p. ( q 


2/3 


3 ) 


„5/3 


( 2 . 20 ) 


see Eqs. (2.6) and (2.17). 


In addition to the solitons, the KdV equation (2.1) 


has delocalized solutions describing periodic waves, the 
cnoidal wavesP The energy (2.2) and momentum (2.4) 


associated with such solutions diverge in the limit of infi¬ 
nite system size unless the waves have vanishingly small 
amplitude. Accordingly, periodic wave solutions of in¬ 
terest here correspond to the harmonic regime when the 
nonlinear term in Eq. (2.1) can be neglected. This leads 
to the relation Cl = — Q 6 between the dimensionless fre¬ 
quency Cl and the wave number Q. In classical mechan¬ 
ics the energy and momentum of such a wave are pro¬ 
portional to the square of its amplitude. In quantum 
mechanics, however, the periodic wave solution corre¬ 
sponds to a phonon with well defined energy and mo¬ 
mentum. Restoring the units, we obtain £xdv = frCl/t* 
and p = hQ/a *. The wave dispersion relation can now 
be converted to the phonon spectrum. It has the form of 


Eq. (2.19) with e(s) given by 


^-phonon (^) ^ 


( 2 . 21 ) 


Use of the classical equation of motion for the evalua¬ 
tion of the excitation spectrum relies on the assumption 
that the quantum uncertai nty o f the field F is negligi¬ 
ble. With the help of Eqs. (2.7) and ( |2.15 1, the magni¬ 
tude of quantum fluctuations of F with the length scale 
~ 1 /A relevant for a classical soliton [see Eq. (2.5)] 
can be estimated as SF ~ 1/A^ ~ A. The condition of 
applicability of the classical description SF <C (F) ~ A 2 
then leads tod> 1, which translat es to Vo 3> 1 for the 
dimens ionles s classical momen tum (2.4) and to p p* 
in Eq. (2.19) or s 1 in Eqs. (2.20) and (2.21). 


For excitations with small momenta p < p* the 
quantum fluctuations can no longer be neglected, and 
Eqs. (2.20) and (2.21) are inapplicable. Fortunately, at 
p —> 0, i .e., d eep in the quantum regime, the KdV Hamil¬ 
tonian (2.10) allows further simplification. Indeed, the 


second term on the right-hand side of (2.10) has a higher 


scaling dimension than the first one, hence its effect on 
long-wavelength excitations reduces to being merely a 
small perturbation.Negle cting th is term and employ¬ 
ing the well-known mapping2iI32I33] between bosons and 
fermions, 


we arrive 


at the fixed-point HamiltoniarfSl 


( 2 . 22 ) 


H le 


2m 


-[ dx\(d x *?{d x *): t (2.23) 

* Jo 


where the symbols • denote the normal ordering with re¬ 
spect to the fermionic vacuum in which all single-particle 
states with positiv e (negative) wave numbers are empty 
(occupied) EH 32 M Application of thejdentity (2.22) also 

yields the fermionic representatiorP^^EH 0 f the momen¬ 
tum operator (2.9), 


P=-ihf dx \ ^(x) d x \ (2.24) 

Jo 

and the relatioip2Ml 


d x § = 27t ■ 4'^(a;)4'(a;) ■ 


(2.25) 


which shows that the held F introduced in Eq. (2.15) 
above is proportional to the density of the effective 
fermions. 


The boundary condition \2 .\\ ) and the relation (2.25) 


imply that SN = f dx\ 'F 1 '(:r)'I , (:r) : = 0. Any eigenstate 
of SN with SN = 0 can be viewed as a superposition of 
the particle- and hole-type elementary excitations. The 
particle excitation is obtained by promoting a fermion 
from the Fermi level (which corresponds to the single¬ 
particle state with wave number zero) to one of the un¬ 
occupied single-particle states, whereas in the hole ex¬ 
citation a fermion is removed from one of the occupied 
single-particle states and placed at the Fermi level. It 
is easy to see that such particle and hole excitations are 
eigenstates of F/fermion and P. The corresponding eigen¬ 
values £fermion and p obviously satisfy 


£ fermion (p) A 


P 

2m* ’ 


(2.26) 


where the +(—) sign corresponds to the particle (hole) 
excitation. This expression can be cast in the form (2.19) 
with 


^fermion(s) — As 


(2.27) 


Note that Eqs. ( |2.20 ) 
|e(s)| ~ 1 at s ~ 1. 


\2.21 ), and ( |2.27 ) yield 
This observation indicates that 
the classical-to-quantum crossover at p ~ p* is the only 
crossover that takes place in the system. 

In view of Eqs. (|2. 15h , (I2.16I) and ( 2.251) , the classical 


soliton excitation corresponds to a hump in the fermionic 
density made up of N ~ No ~ fermions drawn from 
a uniform background. According to the discussion 
above, the classical description is applicable as long as 
N 1. On the contrary, the quantum counterpart of 
the classical soliton, the particle excitation, has exactly 
one excited fermion, which can be interpreted as N = 1. 
Quite naturally, the classical-to-quantum crossover oc¬ 
curs at N ~ 1, i.e., when the discreteness of N can no 

longer be ignored. _ 

In the above consideration we treated a* in Eq. (2.10) 


as a positive parameter. In fact, it can have either sign. 
However, it is easy to see that the operators iFi<dv( a *) 
and iFKdv(~ a*) are related by the transformation 
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<3> —» — $, which amounts to the particle-hole transfor¬ 
mation 'F —» 4/1 for fermions, see Eq. (2.22). Under such 
particle-hole transformation -ffKdv(~ ctTJ— > — ^Kdv( a ») 
and P —>• P. Accordingly, the spectra of the ele¬ 
mentary excitations of -ffKdv( — a*) are related to those 
of i?Kdv(o*) and can be written in the form similar 
to Eq. d2A9k as 


£Kdv(p) = — e*e(p/p*). 


(2.28) 


Note that the particle-hole transformation changes the 
sign of the fermionic density (2.25). Thus, whereas 
the solitonic excitation of -f/Kdv( a *) carries Af > 0 
fermions and corresponds to a hump in the fermionic den¬ 
sity, its particle-hole-transformed twin has Af < 0, which 
amounts to a depression. In nonlinear optics literature 
these two kinds of solitons are often referred to as bright 
and dark solitons, respectively.^ 


III. MODELS 


Instead of attempting to evaluate the excitation spec¬ 
trum of the quantum KdV model (2.10) directly, we 


rely on well-known solvable models, viz. the Lieb-Liniger 
model and the quantum Toda model. In this section, 
we demonstrate that in judiciously chosen scaling limits 
the low-energy theorie s des c ribin g these models reduce to 
that defined by Eqs. (2.7)-(2.10 1 . The reduction hinges 


on the smallness of the parameter 


c = 


p* 


vp* 2 TO*U ’ 


(3.1) 


which characterizes the relative magnitude of the 
KdV cont ribution to the excitation spectra [see 
Eq. (2.12)] evaluated near the classical-to-quantum 
crossover at Although the low-energy Hamilto¬ 

nians for both the Lieb-Liniger and the quantum Toda 
models contain non-universal contributions absent in 


Eq. (2.8), below we show that these contributions affect 


the spectra only in the second order in £, and 

e(p) = vp*[s±(e(s) + 0(( 2 )], s = p/p *. (3.2) 

Here the +(—) sign corresponds to the quant um Toda 
(Lieb-Liniger) model, and e(s) [see Eq. (2.19)] are the 
dimensionless crossover functions describing the spectra 
of the quantum KdV model (2.10). 


A. Lieb-Liniger model 

The Lieb-Liniger inode :" 1 1 •■‘CH] 



describes bosons with contact interaction. We are inter¬ 
ested in the thermodynamic limit when both the number 


of particles N and the system size L are taken to infinity 
with their ratio, the mean density no = N/L , kept fixed. 
The interaction strength is characterized by the dimen¬ 
sionless parameter 1 ^ 7 = c/no- In the weak repulsion 
regime considered here, 0 < 7 C 1. The Hamiltonian 
(3.3) can also be written in the second-quantized form 


H = 


2m 


dx 


(d x ipy (d x ip) + cn 2 (x) 


(3.4) 


where n(x) = (x)ip(x) is the density operator. Here¬ 
inafter, products of quantum fields at the same spatial 


point, such as n 2 (x) in Eq. (3.4), are to be understood 


as being normal-ordered with respect to the appropriate 
vacuum states, cf. Sec. [TT] 

We are interested in excitations with wavelengths much 
longer than the mean interpar ticle distance 1/no- Follow¬ 
ing the standard prescription j 35 l 3f » l we wr jte 


ip(x) = \Jn{x) e 


(3.5) 


Here n(x) is the coarse-grained (averaged over a region 
much larger than 1 /no) particle density, which can be 
regarde d as a continuous function of x, and the field d 

satisfie^ESI 


[n(x),d(y)] = -iS(x-y). (3.6) 


We will assume that d obeys the periodic boundary con¬ 
dition d(x + L) = d(x). This assumption amountto 
restricting one’s attention to excitations near the zero- 
momentum ground state. 


Substitution of ip(x) in the form (3.5) into Eq. (3.4) 
yields 


H = 


2 TO 


dx 


i{dxdf 


1 

4n 


(d x n) 2 


(3.7) 


For low-energy excitations, deviations of the density n{x) 
from its mean value no are small. It is therefore conve¬ 
nient to write n(x) as 

n{x) = n 0 + -d x ip, (3.8) 

7r 

where the new bosonic field p satisfies 

[d x <p,d(y)\ = -inS(x - y) (3.9) 


and obeys the periodic boundary condition. Successive 
approximations for the low-energy Hamiltonian are ob¬ 
tained by substituting Eq. (3.8) into Eq. (3.7) and ex¬ 
panding in powers of d x ip. 

The leading contribution to the low-energy Hamilto¬ 
nian contains operators of scaling di mensio n two and it 
has the standard Luttinger liquid forrr J 33 * 35 -JEZl 

H o = dx[K(d x d) 2 + K-\d x <pf]. (3.10) 


Here 


Hn 0 irhn 0 

v = -V 7 , K = - 

m mv 


7r 

Vi 


= 4 = ( 3 - 11 ) 
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are the sound velocity and the Luttinger liquid parameter 
evaluated in the leading order in 7 <C 1. Instead of ip 
and ■$, it is convenient to introduce the right- and left- 
moving fields 

<P±(x) = (3.12) 

which satisfy \d x ip+,ip-(y)\ = [<p+(a:), d y ip-] = 0 and 
\ip±{x),tp±(y)\ = ±i.Trsgn(x- y). (3.13) 


Writing the Luttinger liquid Hamiltonian (3.101 in terms 


of (fi± reveals its chiral (i.e., diagonal in the basis of the 
right- and left-movers) structure, 


H n = 


hv 

47T 


dx ^2( d x<Pv) 2 - 


(3.14) 


i/ = ± 


In addition to the Hamiltonian, we need the momen¬ 
tum operator P = —ihj'dxip^dx'il). Subs tituting here 
Eq. (3.5) and taking into account Eq. (3.8), we obtain 


P = — hjdxndx'd = — — Jdx {d x tp)(d x d). (3.15) 

Written in terms of the right- and left-movers ip±, see 
Eq. (3.12), the momentum (3.15) takes the form 


P = P+ + P_, P± =±— dx(d x p±) 2 . (3.16) 


Comparison of Eqs. (3.14) and (3.16) shows that the Lut 


tinger liquid Hamiltonian (3.14) can be written as 


H 0 = v(P+ - P_). 


(3.17) 


The nonlinear corrections to the excitation spectra 
arise due to higher-order terms in the gradient expan¬ 
sion of the Hamiltonian (3.7). Collecting the chiral terms 


with scaling dimensions three and four in this expansion, 
we write the resulting contribution in the form of the 
quantum KdV Hamiltonian (2.10), 


H, = 


h 2 


127T to* 


dx [(d x ^) 3 + a*(dlip v ) 2 , (3.18) 


where the effective mass m* and the emergent length 
scale a* are defined by 


to 

to * 


3 

4 VK 1 


a * tiq — 


p- 3/2 

2-7T 2 


(3.19) 


and satisfy to/to* <C 1, a*roo 1. Note that the 
result for the effective mass agrees with the general 


expressio: 


to 1 d(vno) 

to* 2 vVK dn 0 

valid for any Galilean-invariant system. 


(3.20) 


Although the second term in Eq. (3.18) has a higher 


scaling dimension than the first one, its effect on the ex¬ 
citations with wavelengths of order a* [or , equ ivalently, 
with momenta of order p* ~ fi/a*, see Eq. ( 2.18 1 [ is com¬ 
parable with that of the first term in Eq. (3.18). The key 


observation is that for weakly interacting bosons 

37T 


p* 

hriQ 


K 3/2 


< 1, 


(3.21) 


i.e., such excitations belong to the realm of the long- 
wavelength description based on the gradient expan¬ 
sion. Changing the integration variable in Eqs. (3.14) 
and (3.18) to ( = x/a 
H = Hq -f- -)-... as 


we write the expansion 


H — up* (Hq + (hi + + C 2 ^2 + •■■), 


(3.22) 


with the parameter ( introduced in Eq. Using 

Eqs. (13.11 ), (13.191) , and (13.211) , we find 


*= 8 K < L 


(3.23) 


The operators h 0 and h\ in Eq. (3.22) are given by 


(3.24a) 

^ E [(^^) 3 + ( 9 I ^) 2 ], ( 3 - 24b ) 


cf. Eqs. (3.14) and (3.18). The third term in Eq. (3.22) 


includes non-chiral contributions with scaling dimensions 


three and four omitted in Eq. (|3.24b|), (djp±) 2 (d^ip T ) 


and (9?p + )(9?p_). Contributions in Eq. (3.22) that are 


higher order in £ originate in the expansion of the sec¬ 
ond term in Eq. (3.7), the so-called quantum pressure,^ 
in d x <p. The first term in this series yields /12 in 
Eq. (3.22), which contains operators of scaling dimension 
five \d£ip v ) 2 (dzip„). 


The expansion (3.22) allows us to classify various terms 


in the low-energy Hamiltonian according to the order of 
magnitude of their contributions to the energy of chiral 
excitations. It should be emphasized that this classifica¬ 
tion is different from the usual notion of scaling dimen¬ 
sion. Indeed, the latter is relevant for the description 
of excitations in the limit of infinitely long wavelengths, 
whereas we are interested in excitations characterized by 
long but finite wavelengths of order a*. 

We proceed by singling out the first two terms in the 
expansion (3.22), H = vp*(ho + (hi), treating the re¬ 


mainder of the expansion as a perturbation. Taking into 
account that H commutes with P±, we consider a si¬ 
multaneous eigenstate of the operators P+, P_, and H 
with eigenvalues p+ = p ~ p*, p_ = 0 , and e, respec¬ 
tively. In the right-moving sector (i.e., when acting on 
states with p_ = 0 , such as the one we discuss), the un¬ 
perturbed Hamiltonian H coincides with that defined by 
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Eqs. ( 2.7| )-( 2.10), with a* in ffKdv [s ee Eg . (2.10 )] re - 
placed with — a*. Using Eqs. (2.12), (2.28), and (3.1), 
we obtain e{p) = vp * [s — £e(s)] with s = p/p *. 

Perturbation theory in H — H = up* {(h' 1 +( 2 h 2 + ■ ■ ■) 
establishes the correspondence between the eigenstate 
of H we consider and the eigenstate of the full Hamilto¬ 
nian H with the energy e = e+Se. The latter state is also 
an eigenstate of the total momentum P [see Eq. (3.16)] 
with the eigenvalue p. The leading corrections to the 
energy arise in the second order in h\ (note that h\ 
has zero expectation value in the right-moving eigen¬ 
state of H), and in the first order in h 2. Both contri¬ 
butions can be estimated as Se ~ vp^( 2 . Accordingly, 
e(p) = e{p) + 0(up*C 2 )j leading to Eq. (3.2). 


B. Quantum Toda model 

The classical Toda model 21 J describes a chain of N par¬ 
ticles with exponential nearest-neighbor interaction, 


H = H 


+V 0 e- 2D,/a ° 

2m 


Di = xi + 1 - xi . (3.25) 


We consider a system of size L in the thermodynamic 
limit taken at a fixed particle density ng = N/L. The po¬ 
tential energy in Eq. (3.25) is minimized in the ordered 
configuration aq < aq < < aqv, in which the dis¬ 

tance between neighboring particles Di equals its mean 
value 1/no. At finite but sufficiently low energies, the 
deviations from the mean SDi = Di — 1/no remain small 
compared with the interaction range ag . In such a weakly 
anharmonic regime the interaction potential in Eq. (3.25) 


can be expanded in series in SDi/ao . After taking a 
continuum limit and making an appropriate change of 
variables, the leading terms of such expansion can be 
caslPS in the form of the KdV Hamiltonian (2.2). The 


model (3.25) is integrable and supports solitonic exci¬ 


tations with arbitrarily high energies, well beyond the 
weakly anharmonic regime. 21 However, these excitations, 
the Toda solitons, can no longer be described in the KdV 
framework. 

In the quantum Toda modelj- 12 ^® coordinates of the 
particles Xi and their momenta pi are replaced with the 
operators satisfying [xi,pi>] = ihSip. Unlike the clas¬ 
sical Toda model, the very existence of the weakly an¬ 
harmonic regime is not a priori guaranteed even at low 
energies but hinges on the smallness of quantum fluctu¬ 
ations of 5Di. The lower bound SD on the magnitude 
of these fluctuations can be estimated as the amplitude 
of zero-point oscillations of the positions of the particles 
near the respective potential minima. It is convenient to 
characterize the range and the strength of the interaction 
potential in Eq. (3.25) by the dimensionless parameters 


1 


a = 


a 0 n 0 


P = 


1 y/rnV 0 
2 ae“ hnn 


(3.26) 


In terms of these parameters, SD ~ ao/3 -1 / 2 and the 
necessary condition for the existence of the weakly an¬ 
harmonic regime SD -C ag translates to 

/3 > 1, (3.27) 

irrespective of the value of a. _ 

The quantum Toda model is integrableEB !2 0- an( ] its 
properties can be studied analytically-® at arbitrary val¬ 
ues of the parameters a and /3. However, our goal is to 
describe excitations in the universal KdV regime rather 
than to explore various regimes of the quantum Toda 
model. Therefore, we make a simplifying assumption 


a » 1 . 


(3.28) 


The advantage of working in the dilute limit (3.28) is two¬ 
fold. First, as shown below, in this limit ( ~ 1//3, hence 
Eq. (3.27) guarantees the existence of the KdV regime. 
Second, excitations of the Toda model in the dilute limit 
can be studied 12 18 ^by consid ering instead the hyperbolic 
Calogero-Sutherland inodeT^l 1 ® 


d 2 


9m 1 ^ ^ S)x 2 ^ ^ 


A(A-l) 


; a.Q sinh 2 [(xi -xi')/a 0 ] 

(3.29) 

For large A the distance between neighboring particles 
is close to l/?ro ap, s ee Eq. (3.28). Therefore, the 
sinh function in Eq. (3.29) can be approximated by expo¬ 


nential, and the interaction can be restricted to near est 
neighbors. Thus, the model (3.29) is equivalen^^KSI t 0 
the Toda model (3.251) with 


U 0 = 


4fi 2 A(A - 1) 


mas 


Substitution into Eq. ( 3.26| ) gives 

fi = VA(A- 1) e~ a « Xe~ a , 


and the condition Eq. (3.27) yields 


A>e° 


(3.30) 


(3.31) 


(3.32) 


The inequalities (3.28) and (3.32) define the Toda limit 


of the hyperbolic Calogero-Sutherland model (3.29). Un¬ 


like the quantum Toda model (3.25), the hyperbolic 
Calogero-Sutherland model (3.29) is solvable by the 
asymptotic Bethe ansatzp^- which allows it to be han¬ 
dled on equal footing with the Lieb-Liniger model. 

For long-wavelength excitations, pi and Di in 


Eq. (3.251 vary with l on the scale much larger than unity, 


and l can be treated as a continuous variable. Replacing 
the summation over l in Eq. (3.25) with the integration, 


we write the low-energy Hamiltonian as 


H= dl 


P\l) 

2m 


V 0 e 


— 2D(l)/ao 


(3.33) 
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where D(l) = x(l + 1) — x(l), and the fields x(l) and p(l) 
satisfy [x(l),p(l')] = ihS(l — l'). 

The Hamiltonian ( |3.33 ) describes the strongly interact¬ 
ing quantum fluid in terms of the so-called Lagrangian 
variables j 9 l 4 °l 41 l i n which the position of the fluid ele¬ 
ment is specified by l rather than by its physical coordi¬ 
nate x(l). The consideration in Secs. |H| and III A| how¬ 
ever, follows the conventional bosonization schem^I® 7 - 
based on the EuleriarP® 1 description of the quantum 
fluid. A switch to the Eulerian formulation is accom¬ 
plished using 


dl = n(x)dx , 


(3.34) 


where n{x) is the coarse-grained density operator, and 
writing the momentum per particle p(l) as^ 


p(l) = -hd x tf, 


(3.35) 


where the field d has essentially the same meaning as 
that in Eq. (3.51 and satisfies the same commutation re¬ 
lation (3.6). Indeed, substituting Eqs. ( |3.34 ) and ( 3.35 ) 
into the expression P = Jdlp(l), we recover Eq. (3.15) 
for the momentum operator, P = —hf dx n(x)d x d. 

The Eulerian form of the interparticle distance D(l) in 
Eq. (3.33) follows from the identity D(l) = [ e dl — l]a;(Z). 
Replacing di with n~ 1 d x here [see Eq. ( 3.34[>] , we find 


D = (e 
1 


- 1 ) 


X 


11.1 

n 2 \ n x n 


1 U 1 

o! n n n 


(3.36) 


i.e., D depends on the density n( x) and its derivatives. 
Predictably, D reduces to l/n in the long-wavelength 
limit. 


Using Eqs. (3.34| and (3.35), we rewrite the Hamilto¬ 
nian (3.33) as 


H = J dxn(d x d) 2 + U[n], (3.37) 

where the interaction energy U[n\ is a functional of den¬ 
sity, 


U[n] =Vq J dxi 


a —2D/ao 


(3.38) 


with D given by Eq. ( 3.36[ ). Note that the field d enters 
the Hamiltonian (3.37) only via the kinetic energy term, 
which has the same form as that in Eq. (3.7). This is 


a direct consequence of the Galilean invariance of the 


models ( |3.3[ ) and ( |3.25 ). 

The remaining consideration parallels that in 


Sec. Ill A we write the density n(x) in the form (3.8) 


and expand the Hamiltonian (3.37) in d x tp. The leading 


contribution in this expansion is the Luttinger liquid 
Hamiltonian (3.10) with 


K = 


4a 2 P 


< 1 . 


(3.39) 


In the next step, we pick terms containing integrals of 
(dx i p)(d x , d) 2 , ( d 2 (p ) 2 , and (d x p) 3 . Expressin g <p a nd d 
via the right/left-moving fields tp± [see Eq. (3.12)] and 
omitting all non-chiral terms, we obtain the KdV contri¬ 
bution 


H, = 


h 2 


-Jdx [(d x <A;) 3 ~ a*{dlip v ) 2 . (3.40) 


127TTO. 


Except for the sign of the second term on the right-hand 
side, Eq. (3.40) has the same form as Eq. (3.18). The 


effective mass m* and the emergent length scale a* in 


Eq. (3.40) are defined by 


Hi 9 HTt ^ 

— = “7^ = ^. 


(3.41a) 

(3.41b) 


again in agreement with Eq. (3.20). Taking into account 
Eqs. (3.27) and (3.28), we find m/m * 1, a*7io 1. 


Substitution of Eq. (3.41b) into Eq. (2.18) yields 


P* 

Hn-o 




= ^Vk«i. 


(3.42) 


for the crossover momentum. The expansion parameter ( 
[see Eq. (3.1)] is then given by 


c = 


87T/3 27r 2 


a 2 K < 1. 


(3.43) 


Focusing on excitations with wavelengths of order a*, 
we change the integration variable to £ = x/a* and write 
the gradient expansion of low-energy Hamiltonian in the 
form (3.22). The Luttinger liquid contribution in this 


expansion, ho, is given by Ecp (3.24a), and the KdV con¬ 


tribution hi differs from Eq. (3.24b) only in the sign of 
the second term on the right-hand side. The remaining 
terms in the expansion (3.22) are non-universal. The 


third term, h! x , contains non-chiral contributions with 
scaling dimensions three and four, whereas the fourth 
term, h 2 , is composed of the chiral contributions 
(d£ip u ) 2 (d 3 ipv), and (d^(p v ) 2 , which originate in the ex¬ 
pansion of the interaction energy Eq. ( 3.38[ ). Impor¬ 
tantly, various operators enter h\ and h 2 with numer¬ 
ical coefficients that have finite limits of order unity 
at a —> oo. 


Repeating the analysis of Sec. |III A[ we conclude that 
and h 2 contribute to the energy of chiral excitations 
with momenta p ~ p* only in the second order in </. 
Therefore, to first order in / the excitation spectrum is 
given by Eqs. (2.12) and (2.19), leading to Eq. (3.2) with 


a positive sign of the second term on the right-hand side. 

IV. ELEMENTARY EXCITATIONS 

As shown in Secs. [TT] and m exci t ation spectra of 
the chiral model defined by Eqs. (2.7)-(2.10), the Lieb- 


Liniger model (3.3), and the quantum Toda model (3.25) 






















































































are given by Eq. (3.21. The non-trivial parts of the spec¬ 


tra are described by the universal dimensionless crossover 
functions e(s), with asymptotes given by Eqs. (|2.20|) 


and (2.21) in the classical limit s 1, and by Eq. (2.27) 
in the quantum limit s <C 1. In this Section we take 
advantage of the integrability of the Lieb-Liniger and 
the hyperbolic Calogero-Sutherland models. We evaluate 
their spectra to first order in £ and extract the crossover 
functions e(s). 


A. Excitation spectra from Bethe ansatz 


The Lieb-Liniger model (3.3) and the hyperbolic 
Calogero-Sutherland model (3.29) are integrable by 
Bethe ansatz D3HH|: their many-body eigenstates are 
parametrized by sets of N rapidities gi, g2, ■ .., gjv, which 
are similar to the wave numbers of free fermions. Lieb’s 
type I and type II excitation^ can be viewed as, re¬ 
spectively, particle- and hole-type excitations of the cor¬ 
responding Fermi seal 11 ' 14 Their momenta and energies 
are given parametrically by 


p(g) = 2n H 


dq'p{q') 


e(q) = 


dq'a{q') 


(4.1) 


where g > go (g < qo) for the type I (type II) ex¬ 
citations, and g 0 is the Fermi rapidity. The function 
p(g) = i d(q — qt) is the density of rapidities in 

the ground state. In the thermodynamic limit, p(q) 
can be viewed as a continuous function of g, normalized 


rg o 


dqp(q) = n 0 , 


1 -qo 


and satisfying the Lieb equation 11 13 

p (^ + 27r/ dq '®'( q -4)p( q ') = 


1 

2t r’ 


(4.2) 


(4.3) 


where 0(g) is the two-particle scattering phase shift (see 
below). The function a(k) in the second equation in (4.1) 


is the derivative of the energy function introduced in 
Ref. [14j. It obeys the Yang-Yang equatioif 11 12 321 

= 5!f. (4.4) 

m 


l rQo 

v(q) + X- / dq'e'(q - q)a(q) 
2lT J-q 0 


The functions p(g) and a(q) a re, res pectively, even and 
odd functions of their argument (- 11 — Their valu es at the 
Fermi rapidity po = p(qo) and cr 0 = cr(go) satisfj^ 11 12 


— = 2-kHv, 
Po 


&oPo — 


h 2 riQ 

2m. 


(4.5) 


Taking into account the relation between the sound 
velocity v and the Luttinger liquid parameter K [see 
Eq. (3.11)], we obtain 


Po = 


Vk 

2tt 


^0 


= hv\[K. 


(4.6) 


The two-particle scattering phase shift 0(g) in 
Eqs. (4.3) and (4.4) is given 


0 (g) = — 2 arctan(g/c) 
for the Lieb-Liniger model (|3.3[), and b} * 12 l 18 l 


0 (g) = 21m 


lnT A + 


ia 0 q 


— In r ( 1 + 


iaoq\ 

2 ) 


(4.7) 


(4.8) 


for the hyperbolic Calogero-Sutherland mode l ( 3.291. A t 
|g| < A/ao (for a and A satisfying Eqs . (3.28) and (3.32), 
this range includes |g| ~ go), Eq. (4.8) simplifies to 


ia 0 q 


0(g) = dog In A — 2Imlnr^l + ^ 


(4.9) 


The approximation (4.9) is adequate for excitations with 
not too high energies, such as e ~ vp *, and it is 
equivalent 1 2 iM. 2C) I to neglecting the difference between the 
hyperbolic Calogero-Sutherland model and the quantum 
Toda model. For brevity, we shall refer to the results 
obtained in the framework of the approximation (4.9) as 
pertaining to the quantum Toda model. 

The dependence on q in both Eqs. ( |4.7[ ) and (4.9) is 
characterized by a single scale g*, with g* = c for the 
Lieb-Liniger model and g* = 2/ao for the hyperbolic 
Calogero-Sutherland model. In the regimes we consider, 
these scales are small compared with the respective val¬ 
ues of the Fermi rapidities, 


- ~ C « 1, 

<7o 


(4.10) 


where £ is the small parameter introduced in Sec. dm see 
Eqs. (1 3.1 1) , (|3.23|), and (|3.43|). Moreover, it follows from 
Eqs. (3Tll), 13)211), 13.391), (|3.42h, and (0)6]) that 


(4.11) 


p* = 6fyo 0 <7* 


for both models considered. Equations ( |4.1[ ) then 
show that |g — go| ~ g* corresponds to p(q) ~ p* and 
e(g) ~ wp». Thus, in order to study the excitation spec¬ 
tra at the classical-to-quantum crossover, it is sufficient 
to find p and er in the vicinity of one of the Fermi rapidi¬ 
ties, say, at g « go- Accordingly, it is convenient to work 
with the “shifted” dimensionless rapidities 


t = 


q - Qo 

q* 


(4.12) 


instead of q. 

In terms of t, the classical-to-quantum crossover 
at p ~ p* corresponds to t ~ 1, which is well within 
the range |t| <C qo/q* ~ 1/C- At such t, the difference 
between the normalized functions 

<>(«) = cW = (4.13) 

Po 00 

is very small, of order (Cl. However, this difference 
cannot be neglected as it is responsible for the nonlinear 
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corrections to the excitation spectra. Indeed, with the for the quantum Toda model, and the relation 


help of Eqs. (4.11)-(4.131, Eqs. (4.1) can be cast in the 
form 


p* 

e± — vp 
up* 


— = [ digit), 


2 t t r l 

vo{t) = dt ' 0 °( t '' ) = 2S ( t ')’ ( 4 - 21 ) 


(4.14a) 


/ 0 


7T 




= I dt[s(t) - e(t)], (4.14b) 


valid for both models. The relation ( |4.21 ) turns out to 
be applicable at all t, not only at |t| 1. This can be 

shown rigorously, see Appendix [B| 

With Eq. (4.21) at hand, the task of evaluating the 


where 0 < r <C l/( and the + (—) signs correspond to 
the type I (type II) excitation. 


B. Crossover functions and their properties 


crossover functions reduces to finding the function Qo{t) 


introduced in Eq. (4.18a). We therefore turn to the Lieb 
equation (4.3), divide it by po [see Eqs. (4.6) and (4.13)], 


change the variable according to Eq. (4.12), and write 
the resulting equation as 


According to the discussion in Secs. [TT] and III the 
excitation spectra at C —► 0 are given by Eq. ( |3.2 ), which 
we write here as 

£±(p) = up* [s + (e±{s) + 0(C 2 )] , s — p/ p*. (4.15) 


dt'Q'{t-t')g{t') 

^ J-2k 0 


1 


2ttpo 


(4.22) 


Comparison of Eq. ( |4.15[ ) with Eqs. ( |4.14| ) then yields the 
dimensionless crossover functions e±(s) in the form 


where kg = qo/q * ~ 1 /£, see Eq. (4.10). 

For the Lieb-Liniger model, the kernel 0'(f) in 
Eq. ( 4.22| ) is given by 


e±(s) = ±E(±s), 


(4.16) 


B'W = r 


t 2 


(4.23) 


where s > 0 and the function E{S) is defined parametri¬ 
cally by 


whereas 1/po cx C, 1 ^- Accordingly, in the limit C 
Eq. (4.22[) turns into a homogeneous equation, 


S(r) = lfdt0o(t), E(t) = ^£dt m {t). (4.17) Qo{t) + ^ j Jt'Q'{t - t') Qo {t') - 0. (4.24) 


Here r may have either sign, and the functions Qo(t) and 
rjo(t) are given by 

Qo{t) = )im g(t) = lirn q{t), (4.18a) 


Augmented with the condition p 0 (0) = 1, Eq. (4.24) 
has a unique solution for the normalized density go{t)- 
This solution is deri ved i n Appendix [C] At large \t\ 
it agrees with Eq. (4.19), as expected. Note that 

VFr 3/2 

— on 


Vit-^eoV) 


hence the integral on 


Vo it) = lim - [q{t) - g{t)], 


(4.18b) the left-hand side of Eq. (4.24) converges. 


with the limits £ —»■ 0 evaluated at fixed t. 

Further insight is provided b y th e asy mpto tic solutions 
of the Bethe ansatz equations (|4.3|) and (4.4) at q satisfy- 


For the quantum Toda model 

O'(t) = 2 In A — 2Re 1 + it ), 


(4.25) 


ing 9* <C |q — (?o| <C Qo: o r, equ ivalen tly, a t t in the range 
1 |t| <C 1/C, see Eqs. (4.10) and (4.12). As discussed 


above, this range corresponds to the classical regime in 
the excitations spectra. Such classical solutions have 
been found in Refs. [13122] for the Lieb-Liniger model 
and in Ref. [l5] (see also Ref. [T2] for a review) for the 
quantum Toda model. These solutions and the approxi¬ 
mations involved are reviewed in Appendix [A] Substitut¬ 
ing the classical solutions into Eqs. (|4.18[), we obtain 


ffo(t) = 


t > 1, 


14f/7T 1 1//2 , — f » 1 


for the Lieb-Liniger model and 

= / (4</A /2 


Qo it) 


t > 1, 

t|- 1/2 , 1 


(4.19) 


(4.20) 


where ’ipiz) = dlnTiz)/dz is the digamma func¬ 
tion. With Eq. (4.20) taken into account, this gives 
0 '(t — oc |t / | -1 / 2 In \\/t’\, and the inte¬ 

gral in the analog of Eq. ( |4.24 1 would be divergent. 
This difficulty is circumvented by first differentiating 
Eq. (4.22) with respect to t, and only then taking the 
limit C —t 0. The resulting integro-differential equation 

Qo(t) + f dt’Q"{t - t')eoit') = 0 (4.26) 

and the condition po(0) = 1 d efine gpjt) uniquely. 

Although the kernels (4.231 and (4.25) have quite dif¬ 
ferent appearance, their Fourier transforms turn out to 
be cl osely relate d, as are the functions Qo{t) defined by 
Eqs. (4.24) and (4.26). Indeed, as shown in Appendix |D| 


Qo{t) 


Lieb-Linige: 


r = Qo(-t) 


quantum Toda' 


(4.27) 
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Qo(t) 



FIG. 1: The function Qo(t) for the Lieb-Liniger mode l. Th e 
solid line is a plot of the exact result given by Eqs. ( |4.30 1. 
The dotted lines represent Eq. ( 4. 19[ > , applicable at |t| 3> 1, 
i.e., in the classical regime. The dashed lines correspond 
to Eqs. (4.311, which include the leading quantum corrections. 


In view of Eqs. (4.17) and (4.21), Eq. (4.27) translates to 
the relation 


| Lieb-Liniger S) 


quantum Toda’ 


(4.28) 


and Eq. (4.16) yields 


e± ( S ) Lieb-Liniger “ e T( S ) 


quantum Toda’ 


(4.29) 


in agreement with the results of Secs, hh and III 


C. Evaluation of the crossover functions 


Equation (4.24| is of Wiener-Hopf type and can be 
solved analytically, see Appendix 0 The solution satis¬ 
fying po(0) = 1 can be written as 

g 0 (t) = -L= f^sinMr^e-^- 1 ^ 
7TV27T Jo Z 

(4.30a) 

at t > 0 and 


Qo(t) = 


f 


dz 


7rv / 27r Jo z 3 / 2 


1 - 


7TC 


z(ln z—l-\-27vt) 


tan(7T3)r(z) 


(4.30b) 


at t < 0. Simple poles at integer z in the integrand of 
Eq. (4.30b) are understood in the Cauchy principal value 
sense. The function g 0 (t) is plotted in Fig. [T] 

As expected for a solution of an integral equation with 
a non-singular kernel [see Eqs. (4.23) and ( |4.24 )], £>o(i) 
is an analytic function. It is easy to check that Qo(t) > 0 
and dgo(t)/dt < 0. The first of these in equali ties i mplies 
that the function E(S) defined by Eqs. flfTFI ) and ( fL2lj ) 
is analytic. This function is plotted in Fig. [2 It satisfies 
E"(S ) > 0 at all S and E( 0) = E'( 0) = 0. Accordingly, 
the crossover functions e±(s) [see Eq. ( 4.16| >] and their 
derivatives e(j_(s) vanish at s = 0, whereas e' + (s) > 0 and 
e'_(s) < 0 at finite s. 


E(S) 



FIG. 2: The function E(S) for the Lieb-Liniger model. The 
solid line is a plot of the exact result obtained by substituting 
Eqs. (4.30) into Eqs. ( 4. 17| ) and ( |4.21 1 . The dashed line is a 
plot of E (S) = S 2 , applicable in the quantum regime |S| <g 1 
[see Eq. (4.361]. 


The asymptotes of the crossover functions in the clas¬ 
sical (s 1) and quantum (s - Cl) r egimes can be found 
analytically. At large |t|, Eqs. (4.30) yield 


Qo(t )| 
eo(*)| 


L»i = 

1/2 

1---1 

+ 

t-H 

oo 

(5” 

+ 

T—1 

-*»i 

4£ 

7r 

1/2 

' In |87 t<| + 1 

|47rt| . 


, (4.31a) 
, (4.31b) 


in agreement with Eq. ( 4.19| ). Substitution into 
Eqs. (4.17) and (4.21) then gives 


E(S) 


ls»i =^+3 5 - 


(4.32a) 


= g(¥) 2/3|5|5/3 ^^ |5|+ ---- (4 - 32b) 


Interestingly, the logarithmic terms, dominating the cor¬ 
rections to Qo(t) in Eqs. (4.31), do not contribute to the 

expansions (4.32). _ _ 

Taking into account Eqs. (4.16) and (4.29), we obtain 
the asymptotes of the crossover functions in the classical 
regime s> 1, 


2 

e_i_(s) = s + — s + ... , (4.33a) 

<M») = -5(f) !/3s5/3 + S <,+ - (4!!3b) 

for the Lieb-Liniger model, and 

«+(») = 5(t) 2 / 3 »' ,3 -5»+--. ( 4 - 34 ») 

e_(s) = — s 3 — J- s + ... (4.34b) 

O 


for the quantum Toda model. The first terms on the 
right-hand sides of Eqs. (4.33) and (4.34| can be deduced 
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from the solutions of the classical equation of motion, see 
Eqs. (2.20) and ( |2.21 ). The second terms in Eqs. (4.33) 
and (4.34) represent the leading quantum corrections. 

At small |t|, the function Qo(t) can be expanded in 
Taylor series. The first two terms of the expansion read 


Qo(t) 


l‘l«i 


= 1 — 7rf/6 + ... 


(4.35) 


(see Appendix [C| , and Eqs. (4.17) and (4.21) result in 

1 


E(S) 


| S |«1 


= S z 


(4.36) 


which yields 


e±(s) = ±s 2 + -s 3 
for the Lieb-Liniger model and 


e±(s) = ±s 2 - ^s 3 


(4.37) 


(4.38) 


for the quantum Toda model in the quantum regime 

s < 1. _ The leading contributions in Eqs. (4.37) 

and (4.38) correspond to the particle and hole excita¬ 
tions of a gas of free fermions with a quadratic dispersion 
relation, see Eq. (2.26). 


V. DISCUSSION 


of one another [see Eqs. (4.15) and (4.16)], with the clas¬ 
sical expressions (2.20) and (2.21) serving as high-energy 
asymptotes. 


At the lowest energies, the classical treatment in¬ 
evitably breaks down. Instead, the system is best de¬ 
scribed in terms of wea kly interacting fermions with a 
quadratic spectruml^UU These effective fermions can be 
used to characterize the elementary excitations in the 
classical regime as well. Indeed, as discussed in Sec. [TTJ 
the bright (dark) solitons carry a macroscopic number 
of fermionic particles (holes). With the decrease of the 
soliton’s energy, this number decreases, becoming unity 
when the bright (dark) soliton reaches its ultimate quan¬ 
tum fate of being demoted to a single particle (hole) ex¬ 
citation of the effective Fermi gas. 


In this paper we considered in detail the integrable 
members of the two families of systems mentioned above, 
namely, the Lieb-Liniger model (3.3) and the quan¬ 


tum Toda model (3.25). The key advantage of working 


with integrable models is that their excitation spectra 
can be found analytically, yielding explicit expressions 
valid throughout the classical-to-quantum crossover, see 
Sec. |IV| The spectra we found satisfy Eqs. (4.15) 


and (I4.29L as expected for models with quantum KdV- 


type low-energy behavior (see Sec. III). 


Our results can be reformulated in terms of the ex¬ 
act spectra of the elementary excitations of the quantum 
KdV model defined by Eqs. (12.71), 12.91), and 12.101) , 


The regime in which the interaction between parti¬ 
cles dominates the energy of a quantum system is of¬ 
ten referred to as classical. Indeed, many properties of 
a quantum system in this regime can be understood by 
solving the corresponding classical equations of motion. 
Paradigmatic examples of one-dimensional systems ex¬ 
hibiting such semi-classical behavior are identical bosons 
with a weak short-range repulsion and identical parti¬ 
cles, either bosons or fermions, with a strong long-range 
repulsion. Despite the obvious difference between these 
two families of systems, it can be shown that their low- 
energy excitations admit a universal description in terms 
of the quantized version of the so-called first Hamilto¬ 
nian structure 2 -^ of the KdV equation. This observation 
implies, in particular, that the spectra of elementary ex¬ 
citations for the members of the two families are related 
to that of the quantum KdV model (2.10), and, therefore, 
to each other, see Eq. (3.2). 

The relevant classical equation of motion is the cele¬ 
brated KdV equation (2.1). This equation has two phys¬ 
ically and mathematically distinct types of solutions: the 
delocalized periodic waves, corresponding to phonons in 
quantum systems, and the solitons. The classical solu¬ 
tions translate to power-law excitation spectra charac¬ 
terized by different exponents for phonons and solitons, 
see Eqs. (2.19)-(2.21). Importantly, in quantum systems 
these two excitation branches are not independent, and 
their spectra can be viewed as analytical continuations 


£±(p) = ±£»E(=fp/p*)- (5.1) 


Here e± a nd p are eigenvalu es of the operators 7?Kdv [see 
Eq. (2.10)] and P [see Eq. (2.9)], and e* and p* are the 
energy and momentum scales expressed via the param¬ 
eters to* and a* of HKdv accordin g to Eq. (2.18). The 
dimensionless function E(S) in Eq. (5.1) is defined para¬ 
metrically by Eqs. (4.17), (4.211, and (4.30). This func¬ 


tion is plotted in Fig. 2 and its asymptotes are given 
in Eqs. (4.32 ) and (4.36). Because the quantum KdV 
Hamiltonian ( 2.10| is chiral, the two excitation branches 
e± (p) coincide with the bounds on the excitation energy e 
at a given momentum p: e_(p) < e < e+(p). 


Finally, we note that the derivation of the low-energy 
quantum KdV theories for the Lieb-Liniger and the quan¬ 
tum Toda models in Sec. m does not rely on their 
integrability and can be readily adapted for generic 
(i.e., non-integrable) systemsP In such systems £+(p) 
excitation mode no longer represents the exact eigen¬ 
state and, th erefore, has a finite lifetime even at zero 
temperature! 5 ! 9 ! 4 ^!] However, the resulting inelastic 
broadening is parametrically small near the classical 
limit .ISMS] Therefore, we expect our results for the ex¬ 
citation spectra to be applicable to all quantum one- 
dimensional systems with a classical limit described by 
the KdV equation. 





































































12 


Acknowledgments 

We benefited from discussions with A. G. Abanov. 
This work was supported by the US Department of En¬ 
ergy, Office of Science, Materials Sciences and Engineer¬ 
ing Division. We are grateful to the Aspen Center for 
Physics (NSF Grant No. PHYS-1066293) for hospitality. 


Appendix A: Classical regime from Bethe ansatz 


In this Appendix, we review asymptotic 


solutions 12 13 118122] 0 f Bethe ansatz equations (4.3) 


and (4.4) at q*/qo ~ C approaching zero. Here q* = c for 


the Lieb-Liniger model and <7* = 2/ao for the hyperbolic 
Calogero-Sutherland model in the Toda limit. 

It is convenient to work with dimensionless rapidities 
fc = q/q *. After such rescalin g, t he Lieb equation (4.3) 


and the Yang-Yang equation (4.4) become 


If 0 , ) 

p( k ) + ^~ dk'Q (k — k')p(k') = —, (Ala) 

2 ^J-ko 2n 

1 f k ° h 2 n 

a{k) + — / dk'Q'Ck - k')a(k') = — fc. (Alb) 

2 W-fc 0 rn 

The dimensionless Fermi rapidity fco = qo/q* is to be 
determined self-consistently from the normalization con¬ 
dition 


r-fco 


n 0 


dk p(k ) = 

-fco d* 


(A2) 


see Eq. (4.2). In the regimes we consider fco ~ 1/C ^ 1- 


We are interested in the behavior of the functions p(k) 
and cr(fc) at k in the range ||fc| — fco| ^ 1. The momenta 
p(k) corresponding to such k satisfy p{k) p*, i.e., they 
belong to the classical regime in the excitation spectra. 


1. Lieb-Liniger model 

For ko 1 the integrals on the left-hand sides of 
Eqs. ([AT]) are dominated by k! satisfying k 0 — \k'\ 1. 
Moreover, it is natural to assume that both at fco — \k\ 1 
and at |fc| — ko 1 the functions p{k) and cr(fc) vary with 
k on the scale of order ko as no other scale is available. 
Therefore, in order to find p{k) and er( k ) at ||fc| — k 0 \ 1, 

it is sufficient to replace the two-particle scattering phase 
shift 0(fc) = — 2arctan k [see Eq. ( |4.7| )] in Eqs. (Al) by 
its asymptote 


with the pole in the integrand treated as the Cauchy 
principal value. 

Unlike the exact phase shift 0(fc), the approximate 
phase shift 0(fc) has a singularity at k = 0. This sin¬ 
gularity translates to non-analyticities at k —» ± fco in 


the solutions of the approximate Lieb equation (A4). In 


fact, the interior of the interval \k\ < ko does not enter 
Eq. (A4) on an equal footing with its exterior: Eq. (A4) 


can be viewed both as an equation for p(fc) in the interior, 
and as a prescription for extending the solution from the 
interior to the exterior. 

Another difficulty that stems from the singular behav¬ 


ior of 0(fc) is that Eq. (A4) does not determine p(k) at 
|fc| < ko uniquely: if p(k) is a solution, then p + 5p with 
8p oc (k% — fc 2 ) -1 / 2 is a solution as well. The ambigu¬ 
ity can be removed by imposing appropriate boundary 
conditions. For ( -A 0, such conditions read 13 


p{k) 


k = ±(/cq —0) 


= 0. 


(AS) 


The conditions ( |A5| ) are justified by the observation that 
the exact density p(k) near the boundaries of the in¬ 
terval |fc| < fc 0 is parametrically smaller than that in 
its interior. Indeed, with the help of Eq. (A4), p(k) at 
fco — |fc| 1 is estimated as p{k) ~ p{ 0) ~ ko- Substi¬ 


tuting this estimate into Eq. (A2), we find fco ~ 7 3 / 2 


On the other hand, at fco — |fc| < 1 [recall that such fc ar e 
not handled properly by the approximate equation (A4)] 
we have p(k) ~ po, where po ~ AT 1 / 2 ~ 7 -1 / 4 is the exact 
value of p(k) at the Fermi rapidity, see Eqs. ( |4.6[ ) and 
(3.11). Thus, the ratio of p(k) near the boundaries to 
that in the interior is indeed small, of order 7 1 / 4 ~ (' l t 2 
[see Eqs. (3.11) and (3.23)]. 

The solution of Eq. JA4 1 subject to the conditions 
is unique and reads 13 -®! 


1 


(fc 2 -fc 2 ) 


1/2 


|fc| < fci 




fc( fc ) = Y T X \ A 

sgn(k) — (k 2 - fc§) 1/2 , |fc| > k 0 - 


(A6) 


count Eqs. (3.11) and (3.23), we obtain 


Substituting Eq. (A61 into Eq. (A2) and taking into ac- 


2 _2K 

0 s/l * 


9 

47rC 


(A7) 


for the dimensionless Fermi rapidity, in agreement with 
the above estimate. 

We now focus on the realm of the long-wavelength exci¬ 
tations, where the universal KdV description is expected 
to be applicable. For such excitations fc is close to ko. 


0(fc) = 0(fc)L. , = -Trsgn(fc) + - . (A3) and Eq. (A6| can be expanded in powers of [fc - fc 0 |/fc 0 , 

|fc|>>1 * leading to 


With this approximation, Eq. (Ala) simplifies to 


8(i*i 


df AU = A, (A4) 

, k-k' 2 tt ’ 1 ’ 


m = - 

7r 


fco-fc \ 1/2 l/fc 0 ~fc x 3/2 

2fcp ) 2 y 2fco 


(A8a) 
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at k < kp and 


p(k) = 5F 


k — fee 

2kn 


- 1/2 


s fk-ko 


2k, 


o 


1/2 


(A8b) 


at k > kp. Changing the variable to t = k — kp [se e 
Eq. ( 4.12 )1 and taking int o ac count Eqs. (4.6), (4.13), 
and (A7), we rewrite Eqs. (A8) in a more compact form 
as 


g(t) — f dt'Qp(t') + .... 

Po 6 J 0 


(A9) 


where 


gp(t) = 0{— i) 14t/7r| 1 2 + d(t)(Trt) 1 / 2 . (A10) 

Equation ( |A9[ ) represents the expansion in small 
0 ~ (k — kp)/kp of the dominant at ( —s► 0 and fixed (t 
contribution to the normalized density g. Alternatively, 
Eq. ( |A9| ) can be interpreted as an asymptotic expansion 
in small £ of g(t) evaluated at fixed t, with t-dependent 
expansion coefficients found at \t\ 1, cf . Ref FRU. Al¬ 

though in such an interpretation Eqs. (A9) and (A10) are 
applicable only at large |t|, it is reassuring that at \t\ ~ 1 
they yield the correct estimate Qp(t) ~ 1. 


Analysis of the Yang-Yang equation (Alb) follows the 


same steps. Approximating Q(k) in Eq. (Alb) by 0(/c) 
[see Eq. (A3)] and taking into account Eq. (A71, we ob¬ 
tain the equation 


9(\k\ - kp)cr(k) + —~tt dk! 

n dk J_ k 


r(k') 

— fe 0 ' * - *' 


2 / nhvk 


. (All) 


As with Eq. (A4), Eq. (All) does not define a(k) 
uniquely and must be augmented with appropriate 
boundary conditions at k —>• ±kp. At kp — |fc| /g> 1 we 
have | cr(fc) | ~ Hvkp ~ fiwy -1 / 2 , which for 7 ~ C, 2 1 

is parametrically larger than cr(k) at k approaching k 0 , 
a(k) ~ (To ~ fiuy -1 / 4 , see Eqs. ( |4.6[ ) and (A7). This 
observation leads to the condition 


CT W|fc = ±(fco-0) =0 > 


(A12) 


cf. Eq. d A5l) . T he solution of Eq. (All) satisfying the 
conditions ( |A12 1 is given bjiH 


°{k) = x 




3/2 


|fc| < k\ 


0 ) 


sgn(fc)^(fc 2 - klf /2 , |fc| > kp. 


. (A13) 

Expanding Eq. (A13) in (k — kp)/kp ~ Qt and taking into 
account Eqs. (4.6) and (4.13), we find 


<;(t) — — 1 — — Qp{t) H—— ' - f dt'gp(t') + ... (A14) 

cto 6 J 0 


with Qp(t) given by Eq. (A10). Similar to Eq. (A9|, 
Eq. (|A14|) can be viewed as an asymptotic expansion in 


small C of the function $(t) evaluated at fixed t such that 
\t\ /§> 1, which corresponds to the classical regime in the 
excitatio n spect ra. Substit ution of Eqs. (A9) and (A14) 
into Eq. (|4.18b|) yields Eq. (|4.21|). 


2. Quantum Toda model 

Instead of attacking the quantum Toda model directly, 
we consider the hyperbolic Calogero-Sutherland model 
( 3.29[ ) in the Toda limit defined by Eqs. (3.28) and (3.321. 
As shown in Ref. mm and confirmed below, under 
these conditions the dimensionless Fermi rapidity kp in 
the Bethe ansatz equations (Al) satisfies 


1 < kp < A. (A15) 

Accordingly, for |/c| ~ kp <C A the two-particle phase shift 
@(fc) = 2Im[lnT(A + ik ) — lnT(l + ik)~\ (A16) 


[see Eq. (4.8)] can be replaced with 

0(^)||fc|«A = hr A — 2ImlnT(l + ik). 


(A17) 


This approximation ignores the existence of the regime 
\k\ /g> A, absent in the true quantum Toda modelj^ but 
is adequate for our purpose. On the other hand, similar 
to the Lieb-Liniger model, the dominant contribution to 
the integrals on the left-hand sides of Eqs. (Al) comes 


from k' satisfying kp — \k/\ 1. Thus, in order to find 

p(k) and cr{k) at k in the range ||fc| — fco| 1, the phase 
shift can be further approximated by 


0 (fc) = Q(k)\ 




= 2fc(ln|A/fc| + l). (A18) 


The Lieb equation (| A1 a) then assumes the forn i 12 ! 18 ! 


p(k) + - 

7r 


dk' In 


-k 0 


A 


k — k' 


p(k') = —. (A19) 


Similar to Eq. (A4) for the Lieb-Liniger model, 


Eq. (A19) serves simultaneously as an equation for p(k) 


in the interior of the interval \k\ < kp, and as a prescrip¬ 
tion for extending p{k) from the interior to the exterior. 
Moreover, at |fc| < kp the firs t ter m on the left-hand side 
ofEq. (A19) can be neglectediASESlIndeed, at kp— \k\ 1 

this term is of order p(k) ~ p(0 ), whereas the second term 
is parametrically larger, of order kp ln(A/fco)p(0) p(0). 

With these appr oximations the Lieb equation can be 
solved exactly)A21iSI The solution reads 


P(k) = x 


1 


v-i/2 


and (3.43) then yield 


kp = 2Xe~ a = 


2a 2 K 


3 

47T^ 


1*1 < 

kp, 

1*1 > 

kp. 


(A20) 

Eqs. 

(3.39 


(A21) 
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for the dimensionless Fermi rapidity. Since A > e“ > 1 
[see Eqs. ( 3.28 1 and (3.32)], k 0 indeed satisfies the 
inequalities ( |A15 ). In the same approximation, 
0(fc) « 0(fc) [see Eq. (|A18|)], the Yang-Yang equa¬ 
tion |Alb| yieldfPlASI 

’°' 1/2 |fc|<fc 0 , 


=4 x 


| -sW- 4 ' 

[ sgn (k)(k 2 - k%) l/2 


(A22) 


\k\ > k 0 


Focusing on the long-wavelength excitations, we 
expand Eqs. (A20) and (A22) in powers of small 


(k — k 0 )/k 0 = 47 r£t/ 3 . With Eqs. (4.6) and (4.13) taken 
into account, the first two terms of these expansions can 
be written as 

g(t) = g 0 {t ) dt'g 0 (t') + ..., (A23a) 


7tC 


?(f) = g 0 (t ) + 
with go (t) given by 

eo{t) = 0(-i)kr 1/2 


dt' g 0 (t') + 


8 (t)( 4 t/T 


A /2 


(A23b) 


(A24) 


Similar to their Lieb-Liniger model counterparts (A9) 


and (A14), Eqs. (A23) can be interpreted as asymptotic 


expansions in small ( of the functions g(t) and c(t) evalu¬ 
ated at fixed large t. Although the numerical coefficients 
in front of the second terms on the right-hand sides of 
Eqs. (A23) differ from those in Eqs. (A9) and (A14), the 


difference does not affect the form of the function 770 (t) 
[see Eq. (4.18b)], which is again given by Eq. (4.21). 


Appendix B: Beyond the classical regime 


Equation (4.21) can be obtained by substituting g(t) 


and <;{t) in the form of Eqs. (A9) and (A14) for the Lieb- 
Liniger model a nd Eq s. (A23) for the quantum Toda 
model into Eq. (4.18b). However, such a derivation is 
valid only at |t| 1 , i.e., in the classical regime in the 

excitation spectra, whereas we are interested in |t| ~ 1 , 
which corresponds to the classical-to-quantum crossover. 
It turns out that relaxing the restriction |t| 1 requires 

merely a replacement of the functions go(t) in Eqs. (A91, 
(A14), (A23), and (4.21|) with the exact solutions of 


Eqs. (4.24) or (4.26). In this Appendix, we derive such a 


generalization of the expansion (A9) for the Lieb-Liniger 
model. Generalizations of Eqs. (A14) and (A23) can be 
obtained in a similar manner. 


It is convenient to write the Lieb equation (4.22) as 

1 


£-2k 0 [e] = 


2ttpo 


(Bl) 


Here k 0 ~ 1/C, Po ~ C 1/2 [see Eqs. ( |4.6| ) and ( |A7[ )], and 
the functional C T is defined by 

Cr [/] = m + ^ fdt'Q'it - t')f(t'), (B2) 


where 0(f) = — 2arctant. The solution of Eq. ( |B1[ ) can 
be viewed as a function of two variables, t and (. We 
are interested in the behavior of g(t, C) at arbitrary t and 
small C <C nrinjl, |f| -1 }- 


Replacement of the phase shift 0(f) in Eq. (B2) with 


0(f) given by Eq. (A3) leads to the approximate Lieb 
equation 


^2k 0 [g] — 


1 


27rp 0 


(B3) 


which, unlike Eq. |B1J), is exactly solvable. To avoid con¬ 
fusion, in this Appendix we use a tilde to distinguish 
the e xact solution g(t. C) of the approximate Lieb equa¬ 
tion (B3) from the solution p(t , C) of the exact Lieb equa¬ 
tion (Bl). At C|f| 1, the function g{(,t) can be ex¬ 

panded as 


_ _ __ 

Q(t,0 = 6o(t) + -fTQiit) + 
6 


(B4) 


with gi (t) = f dt'g 0 (t’) [see Eq. (A9)]. The function 
p 0 (f) in Eq. (B4| satisfies the equation 


£oo[f?o] — 0 

and is given by Qo(t) = 0 (— t)| 4 t/ 7 r 
Eq. dMOk]. 


(B5) 

1 ^ 2 + 0(t)(7rf ) -1 / 2 [see 


We seek the solution of Eq . (|B1| in the form of an 
expansion inspired by Eq. 


e(t, 0 = 80 (t) 


4 J 

7rC 

6 


8i(t) 


(B 6 ) 


where the functions 80 (t) and £>i(f) are independent of C- 
Since g( 0 , C) = 1 for any C, these functions satisfy 


eo(Q) = i, f?i(o) = 0. 


(B7) 


On the other hand, we expect the functions g(t, C) and 
g(t, C) to match at 1 <C |t| 1/C in every order in (. 

This leads to the relations 


lim 


80 (t) 


= 1 , 


lim 


|t|->oo 8o(t) ’ Id 

so that g[ (f) = dgi (t) /dt satisfies 

8[{t) 


8i(t) 

8i{t) 


= L 


lim 

\t\ OO 


8o(t) 


= 1 , 


(B 8 ) 


(B9) 


where we used g[(t) = go(t). 


Formally, the right-hand side of Eq. (B 6 ) represents 


an asymptotic expansion of p(t, C) in small C evaluated 
at fixed t. The expansion is applicable as long as the 
second term is small compared with the first one. With 
Eqs. 


and (B 8 ) taken into account, this yields the 
condition C •C min{l, |f:|~ 1 } on C> which translates to 
|t| 1/C for t. Taking advantage of this inequality, we 

introduce the intermediate scale r satisfying 


max {l, |f|} «r« 1/C, 


(BIO) 
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f- 


2ko 


dt'Q'(t — t')g(t') in two, 


/- 2 fcn d< '[- ' '] = I-r dt 'l ■ •] + f-L„ dt 'l ■ •]> and replace 0 


break the integral in 
r ° 2ko dt>[...} = f°_ T ^,., , j_ 2ko , 

and p in the second integral here with 0 and g. (The rel¬ 
ative error of such an approximation is of order 1/r -C 1.) 
With the help of Eqs. (|B1|)-(|B3|), we obtain the equation 


C-Ae] = £t[q]- 


(BH) 


The functions go(t) and g\(t) in Eq. (B6) can now be 
deduced by replacing g and g in Eq. (Bill with the ex¬ 


pansions (B4) and (B61, respectively, and considering the 


behavior of the resulting equation at small C and large r. 
For the equation to hold at large, b ut stil l finite r and 
arbitrarily small ( <C 1/r [see Eq. (BIO)], it must be 
satisfied in every order in ( separately. This observation 
leads to the equations 


£"r [go] — 

£t[pi] = C T [gi]. 


(B12a) 

(B12b) 


Differentiating these equations with respect to r, we re¬ 
cover the relations (B8) and (B9). 

In the limit r —> oo, Eq. (B12a) yields 


£oo[go] = 0, 


(B13) 


where we took into acco unt Eq. (B5). Equation (B13) 
coincides with Eq. ( |4.24 ). The solution of this equation 
satisfying the condition p o (0) = 1 [see Eq. (B7)] is unique 
an d beh aves at l arge |f | as prescribed by the first equation 

in (B8), see Sec. IV C and Appendix C]_ 

Because the integrals in Eq. (B12b) diverge at 


r —^ oo, we first differentiate both sides of this equa¬ 
tion with respect to t and integrate by parts using 
gi(0) = gi (—0 ) = 0. Taking now the limit r —>■ oo and 
using Eq. (|B5|, we obtain the equation 


^oo [pi] = 0. 


(B14) 


Its solution subject to the condition (B9) reads 
g[(t) = go(t). Combining this result with the second 


equation in (B7), we finally arrive at 


Qi{t) = [ dt'g 0 (t'). (B15) 

Jo 


B y con struction, Eq. (B6) wit h gi{t) given by 
Eq. ( |B15 1 has the form of Eq. (A9). We have verified 
that Eqs. (B6) and (B15) with pole approximated by its 


asymptote (4.31b) are in agreement with the expansion 
derived by a different method in Ref. [46] , 


Appendix C: Solution of Eq. (4.24) 


In this Appendix we employ the Wiener-Hopf tech- 
niqu e (see , e.g., Ref. S3) to construct the solution of 
Eq. (14.241. 


Substituting go(t) in the form 
Bo (t) = g+(t) + g-(t), g±(t ) = 9(±t)g 0 (t) 


into Eq. (4.24), we obtain the equation 


with 


/ OO 

dt'Q(t — t')g_{t') = 0 

-OO 


G(t) = 6(t) + %^ = 6{t) - ' 


2tt 


7 r(l -I- t 2 ) ' 


(Cl) 


(C2) 


(C3) 


Upon the Fourier transform, Eq. (C2) assumes the form 
p+(w) + <5(w)p_(u;) = 0. (C4) 

The functions 


/ OO 

dt e lu>t g± (t) 

-OO 


(C5) 


in Eq. (C4| are analytic at ±Imu > 0. Their behavior 


at large |w| is obtained by substituting go(t) in the form 
of the Taylor series gp(f ) = 1 + gp(0)t + ... [recall that 
p 0 (0) = 1] into Eqs. (|C1| and (C5), which yields 


e± ^U|»i ± w 


1 + Po(0)-b • • • 

UJ 


The kernel Q(ui) in Eq. (C4) is given by 

/ OO 

dte iuit g(t) = 1-e-H 

-oo 

It can be factorized as 


GH = - 


F_( W ) 


with 


F+(u) = Q+(u)f+( w), A_(w) = 


g-M 

/-M' 


Here 


g+M = 


p 27r/4 


{uj + ioy / 2 ’ 


g-M = 


0 —3in/4 


(uj — zO ) 3 / 2 


(C6) 


(C7) 


(C8) 


(C9) 


(CIO) 


are Fourier transfor ms o f g±(t) = 9(±t)g 0 (t), where 
Po(i) is th e solution (A10) of the approximate Lieb equa¬ 
tion (B5). The functions g± (to) are analytic in the com¬ 


plex plane with the branch cuts running from iO to 
=F zoo. The functions f±(u>) in Eq. (|C9[) are given by 


/N expjT^Snf 5 ^) - 1 =F ifl } 

/±(w) = 1 2nL / 27r \ - 2JJ , Cll 

r(iT^) ; 

with the same branch cut structure. These functions ap¬ 
proach 1 at lo —> 0, whereas at large cc application of 
the Stirling formula gives the asymptotes 


/±M 


0 ± 27r/4 


M»i 


ITT 


(C12) 
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Because g± (w) and f± (w) have no singularities or z eros 
at ± Imw > 0, the functions F± (w) given by Eq. (|C9[) are 
analytic in these regions. Standard arguments^ 7 ' then 
show that the function 


into Eq ( |4.26 ) and using 0"(f) = — 0"(—t) [see 
Eq. (4.251], we obtain the equation 


Hu) = 


g+M = e-M 

F+(w) F-(w) 


(CIS) 


is analytic in the entire complex plane. Using Eqs. (C6), 
(C9), (CIO), and (C12), we find ■^ r ( w )|| tl ,|_ >00 = 1- There¬ 
fore, T = 1 at all w, and Eq. (C13|) yields 


i r°° 

Qo(t) + ^nj 0 dt'Q"{t-t')e 0 (t')=Q. (D2) 

Writing Qo(t) as 

&o(t) = e+(t) + e±{t) = 0{±t)Qo(t) (D3) 


g±(uj) = F±(uj). 


(C14) 


[cf. Eq. (Cl)], and using the continuity of Qo(t) at t = 0, 
we find 


At small w, Eq. (C14) reduces to 

g±(uj) || w | ^ 1 = Q±(oj), hence at ^f] 1 the inverse 

Fourier transforms 


Here 


, . diO —quit / \ 

Q±{t)= 6±(u) 


(C15) 


eL(t) + dt'g'(t~t')e+(t') = o. (D4) 


G{t) = S(t) - - Ref/'(1 + it), (D5) 


reproduce the classical asymptotes Eq. (4.19). At large w, 


Eq. (C14) yields the expansion (C6) with 


where ^(z) is the digamma function. Fourier transform 
of Eq. (|D5| reads 


2o (0) = - 7r / 1 2 3 4 5 6 ^ 


(C16) 


G{u)= dte lult g(t) = 


1 


1 


1 — e~l w l G(u>) ’ 


(D6) 


resulting in the Taylor series At intermediate w with g^) g i ven by Eq. dC7». With the help of Eq. pil, 


and t no further simplifications are possible. Instead, we 


deform the integration paths in Eq. (C15) to run along 


the respective branch cuts and change the integration 
variables to z = =h iu>/2 n. Taking into account Eq. (Cl), 
we arrive at Eqs. (4.30). 


Appendix D: Derivation of Eq. (4.27) 


Fourier transform of Eq. (D4) can then be written as 

p+(w) + Q{u)q-{w) = 0. (D7) 

Carrying out the inverse Fourier transform, we obtain 

/ OO 

dt'g(t - t')g-{t') = 0, (D8) 

-CO 


In this Appendix, we derive the relation (4.27) between which is the same as Eq. (C2). The corres pond ing equa- 


the functions go (t) for the Lieb-Liniger model and for the 
quantum Toda model. 

Substituting 


0 o (t) = 0o (~t) 


(Dl) 


tion for go(t) then coincides with Eq. (4.24) for the 
Lieb-Liniger model and is subject to the same constraint 
p(0) = 1 at t = 0. Therefore, go(t) for the quantum Toda 
model is identical to go(t) for the Lieb-Liniger model, 
which leads to Eq. (4.27). 
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